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Abstract 

We consider a central system which is coupled via dephasing to an 
open system, i.e. an intermediate system which in turn is coupled to 
another environment. Considering intermediate and far environment 
as one composite system, the coherences in the central system are 
given in the form of fidelity amplitudes for a certain perturbed echo 
dynamics in the composite environment. On the basis of the Born- 
Markov approximation, we derive a master equation for the reduction 
of that dynamics to the intermediate system alone. In distinction to 
an earlier paper [arXiv: 1502.04143 (2015)] where we discussed the 
stabilizing effect of the far environment on the decoherence in the 
central system, we focus here on the possibility to use the measurable 
coherences in the central system for probing the open quantum dy¬ 
namics in the intermediate system. We illustrate our results for the 
case of chaotic dynamics in the near environment, where we compare 
random matrix simulations with our analytical result. 

Keywords Quantum Loschmidt echo, Fidelity, Open quantum system, 
Master equation, Random matrix theory 
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1 Introduction 


Loss of fidelity and decoherence are the twin obstacles to snccessfnl appli¬ 
cations of qnantum information devices. Theoreticians like to consider the 
two separately, while in practical situations both will have destructive ef¬ 
fects on quantum information flow. In this spirit it is important to propose 
a functional dehnition for the hdelity or better the hdelity amplitude of an 
open system. To achieve this we shall build on previous work. Almost 20 
years ago a measurement of the hdelity amplitude in a quantum system was 
proposed [1], though the word itself was not mentioned. A more detailed 
presentation for a kicked rotor followed years later [2], simultaneously with 
an extensive discussion of the theoretical framework [3]. The experiment for 
the kicked rotor was performed successfully in Ref. [4]. The basic idea is that 
pure dephasing between a qubit and some system will cause the off diagonal 
elements or coherences of the density matrix of a pure superposition state 
of the qubit decay like the hdelity amplitude of the remaining system. Our 
proposition is to follow exactly the same reasoning in the case where the 
remaining system is open, and thus interpret the decaying coherences as a 
generalized hdelity amplitude of the intermediate system. 

For this purpose we have to consider a situation discussed in a previous 
paper [5], where a central system coupled to two nested environments is an¬ 
alyzed to determine the ehect of a far environment on the coherence of a 
central system not interacting directly with the far environment. The inter¬ 
mediate system or “near environment” would then be the open system we are 
analyzing using the qubit as a probe and as the source of the perturbation 
causing the hdelity decay in the absence of the far environment. 

Our paper will thus fulhll a double purpose: On the one hand we shall 
introduce a generalized hdelity amplitude and study its behaviour in a “quan¬ 
tum chaotic” setting. On the other we shall deepen the understanding of the 
stabilizing ehect on the central system of the coupling of the near environ¬ 
ment to a far environment. Note that in a previous paper [5] we introduced 
pure dephasing between the central system and the near environment, as a 
simplifying approximation to facilitate calculations, but for the present con¬ 
text this is an essential feature needed in order to establish the relation to 
hdelity decay of a closed system. 

Considering near and far environment as one closed system, the hdelity 
amplitude may be expressed as the expectation value of an echo operator [6], 
which describes the forward and backward evolution of an initial state with 
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somewhat different Hamiltonians. In the present paper we shall derive a 
master equation which describes this echo dynamics reduced to the near 
environment alone. This master equation has the typical structure of master 
equations of Kosakowski-Sudarshan-Lindblad form [7, 8] (henceforth referred 
to as Lindblad equation) and reduce to Lindblad form when the coupling to 
the central system becomes zero. 

We then consider the random matrix model already used in Ref. [5]. In 
that case, the master equation becomes very simple, and allows to obtain a 
closed integral equation for the generalized fidelity amplitude. We compare 
the solution of the integral equation with the generalized fidelity amplitude 
obtained from numerical simulations, and thereby show that the new integral 
equation applies for a very broad range of coupling strengths between near 
and far environment. While the random matrix model used here, is much 
simpler than the one used in Ref. [5], its effect on the fidelity amplitude is 
typically indistinguishable. 

In the next section we shall fix notations and obtain the general master 
equation for describing the reduced echo dynamics in the near environment. 
At the end of this section, we derive the master equation for the random 
matrix model, which is largely equivalent to the model considered in [5]. In 
Sec. 3 we compare the generalized fidelity amplitude obtained from the new 
integral equation with numerical simulations, and in the last section we draw 
conclusions. 


2 Model 

In this section we describe the model in general. The full system consists of 
three parts, the central system, the near environment and a far environment. 
In the following section (Sec. 2.1) we start from a Hamiltonian formulation, 
and perform partial traces in order to obtain the reduced dynamics either 
in the central system or in the near environment. Assuming a dephasing 
coupling between central system and near environment, we can obtain the 
temporal evolution of the reduced density matrix in the central system, in 
terms of an asymmetric unitary evolution (perturbed echo dynamics) in the 
composite environment. In Sec. 2.2 we then use the standard Born-Markov 
approximation to trace out the far environment, and thereby arrive at a 
master equation for the reduced echo dynamics in the near environment 
alone. 
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2.1 Hamiltonian formulation 

For the coupling between central system and near environment, we assume 
that it is of the dephasing type, i.e. that it is given by a single product term, 
where the factor acting in the central system commutes with the Hamiltonian 
describing the dynamics in the central system. Thus, the Hamiltonian for 
this part is given as 

Hc,e = he <S) le + '^c ^ He + Vc <S) Ve , where [he, Uc] = 0 . (1) 

This leads to the following expression for the time evolution of the reduced 
density matrix in the central system [1, 3] 

which yields for its individual matrix elements (coherences) 

= e,\(0) /,(() , /,(i) = . (3) 

Here, we use the common eigenbasis of he and Ve to express Qeit), and there¬ 
fore the energies Sj and Sk are simply the corresponding eigenvalues of he- 
In what follows we will focus on only one such matrix element, and therefore 
suppress the indices j, k from now on. We thus set 

Hx = He + We = He + Vj W and Ho = He + Vk 14 , (4) 

which implies that X = uj — Uk- This shows that under dephasing coupling, 
the decoherence in the central system is given by the decay of the hdelity or 
Lohschmidt echo in the near environment. Turning the argument the other 
way around, this shows that it is possible to measure hdelity amplitudes, 
by coupling the system of interest (i.e. the near environment) to a probe 
system, which at the same time provides the perturbation. Such experiments 
have been proposed and recently realized in different settings using atom 
interferometry [1, 2, 9, 4]. 

Including the far environment We now extend the model to include a 
far environment. We assume the far environment to be as simple as possible 
and that it allows to be taken into account implicitly in the form of a quantum 
master equation. Hence, we write for the Hamiltonian of the full tripartite 
system: 

Hx,r = He,e ® If + Ic ® {le® Hf + -fVe^Vf) , (5) 
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where F = 27r TVg 7 ^/ {hdf). Here, and Vf are normalized in such a way that 
7 ^ gives the magnitude squared of a typical matrix element of the coupling 
term between near and far environment. With df being the average level 
spacing (i.e. the inverse level density) in the spectrum of we see that F 
is just Ng times the corresponding Fermi golden rule transition rate [10, 11]. 
Finally, F itself is related to the decoherence rate for a superposition of Ng 
states in the near environment. In the simplest case (see Sec. 3), 2F is 
precisely the decay rate of the purity in the intermediate system, in the case 
where the coupling to the central system is set to zero. 

Within the Hamiltonian model described by H\^r, the coupling to the 
far environment requires the following modihcation to the expression for the 
hdelity amplitude given in Eq. (3): 

fx,r (t) = tre,f [ Pe,f (0) ‘/'^ ] . (6) 

Differing from the standard formalism, the unitary operators on the left and 
on the right hand side of the initial state are different. This is why the trace 
may decrease in time, leading to the loss of coherence for superposition states 
in the central system [5]. 

2.2 Master equation for the echo dynamics 

We follow the standard derivation of the Born-Markov approximation, e.g. 
Sec. 3.3 of Ref. [12]. However, due to the asymmetric unitary transformation 
implied in Eq. ( 6 ), the following derivation requires some care. Let us denote 
the solution in the Hilbert space of near and far environment as 

and thereby introduce X{t) as the solution in the interaction picture with 
respect to the coupling between near and far environment. From the von 
Neumann equation for Pe,f(t), 


ihdt Pe,f(t) — Hx^r Qe,iit) — Pe,f(t) Ho,r , 


( 8 ) 


we obtain 

ihdtXit) = 7 [1/A(t) X{t)-Xit) Voit)] , Vxit) = Vim . 

(9) 
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The aim of the Born-Markov approximation (to be worked out next) consists 
in obtaining a master equation for the reduced dynamics in the near envi¬ 
ronment alone. That is, an evolution equation for Pe(t) = trf[pe,f(t)] which 
may not be considered a real quantum state, since for the reasons discussed 
above, Pe(t) is neither Hermitian nor trace-preserving. This quasi-density 
operator is related to X (t) as follows: 

Pe(t) = trf X{t) = tTf[Ux{t) < 8 ) Ut{t) X{t) Uo{t)^ 0 , 

( 10 ) 

where U\{t) = and Ui{t) = since H\fl and hfo,o are both 

separable operators. Using the identities 

trf [A (g) 1 X ] = A trf [X ] , trf [X i? 0 1 ] = trf [X ] B 

trf [l 0 yl X ] = trf [X 1 0 A ] , (11) 

we hnd 

Pe(t) = Uxit) Qe{t) Uo(f)^ , Qe{t) =trf[X(f)] . (12) 


Born-Markov approximation Now, we will formally integrate the dif¬ 
ferential equation for X(f), Eq. (9), and plug the result back into its right 
hand side: 


X(t) = e..f(0) - T / dr (y (r) A'(t) - A'(t) ^(r)) , 


h 


ihdtX{t) ='j Vx{t) Pe,f(0) - ^ / dr (UA(r)X(r)-X(r)Uo(r)) 


-7 


2e,f(0)-^/ dr (UA(r)X(r)-X(r)Uo(r)) Vo{t) 


h 


= 7 


14(t) Pe,f(0) — Pe,f(0) Vo{t) ] 


h 


2 rt 


dr [Vx{t) (UA(r)X(r)-X(r)Uo(r)) 


UA(r)X(r)-X(r)Uo(r)) Uo(4 ] 


(13) 


The next step consists in taking the partial trace with respect to the far envi¬ 
ronment. In doing so, we will assume that the first term, which comes before 
the integral over r, will not contribute, i.e. trf [ 14(4 2e,f(0) — Pe,f(0) Vo{t) ] = 
0. While this is true in the random matrix model, in general it might be 
necessary to take this term into account. Even then it does not present any 
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difficulty, since the term is known before hand. With ^e(^) = trf[W(t)], we 
hnd 


= J dr trt[V,(«) ( V,(t)X(t) - Vi(T)) 

- ( n(T) X(t) - X(t) i>„(r)) Ut) ] . (14) 

Here, we perform the two crucial approximations: (i) the Born approxima¬ 
tion, which assumes that the influence of the near environment on the state of 
the far environment is negligible: X{t) ~ ^e('r) ® ^>£(0), and (ii) the Markov 
approximation, that the state of the near environment Pe('r) is changing 
slowly on the time scale of the correlation function Vx{t — s) V\{t). We finally 
assume that the Vx(t — s) Vx(t) quickly approaches zero as s increases, so that 

dt ge{t) = - ^ j dr Hf [Vx{t) ( Vx{t) g^it) ^ ft(0) - ge{t) 0 ft(0) Vo{t) ) 

- ( Vxir) Qeit) (8) ^f(O) - ^e(^) <8) ^>£(0) Vo{t) ) Vo{t) ] 
ds trf [Vx{t) {Vx{t- s) Qeit) <8) ^>£(0) - ^e(^) <8) ^£(0) Vo{t - s )) 

- ( Vx{t - s) Qeit) (g) ^£(0) - Qeit) (g) ^f(o) Vo{t - s) ) Vo{t) ] . 

(15) 

This equation is the equivalent of the so called Redfield equation [13]. Next, 
we will consider each term separately and take advantage of the fact that 
Vx(t) is a tensor product operator with respect to the near and the far envi¬ 
ronment: 

Vxit) = Uxity V:Uxit) (g) Vi U,{t) , 

such that due to trf[H ® B] = A trf[i?]. It is now natural to assume ^?£(0) to 
be diagonal in the eigenbasis of i.e. [^)f(0),hf£] = 0 (such would be the 
case for a thermal state). Then we find 

dt Ut) = -^ I ds Cis) { Uxit)^ V: Ux(s) Vj Ux{t - s) g,{t) 

- Uxity v: Uxit) Ut) Uoit - 5)t Vi U,{t - s) 

- Ux{t - 5)t v: Ux{t - s) ~g,{t) Uo{t)^ Vi Uo{t) 

+ ~ge{t) Uo{t - s)t Vi UxisY Vi Uo{t) } , (16) 
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with the real function C(s) = tif [f/f (s)l 1/f 17f(s) Vf ]. Going back to the 
Schrodinger picture, we obtain 

ihdt Qe{t) = Hx Qe{t) - Qe{t) Hq 

- ^ { K' Ta Qeit) - v: g,{t) To - Ta Pe(t) v: + Pe(t) Tq K' } , (17) 

where 

poo 

Tx= dsC{s)Uxis)V:Uxis)K (18) 

Jo 

In the simplest case, G(s) may be approximated by a delta function, such 
that Fa = CoV^/2, where Cq is the area under the function G(s), and the 
factor 1/2 comes from the fact that the integral goes only over the positive 
half axis. In this case, the master equation (17) is of Lindblad form with the 
Hermitian Lindblad operator \/CoV/. Note however, that there are many 
interesting cases, where C{s) has a different functional behavior, which then 
yields a rather unusual dissipation term. 

RMT model We here use the approach of constructing the average den¬ 
sity matrix proposed in [14] rather then calculating properties for each mem¬ 
ber of the ensemble and averaging these properties afterwards. For details 
see [15, 16]. Let us assume the master equation derived above is of Lindblad 
form with a single Lindblad operator V/. Then, we may choose Cq = 1 with¬ 
out restriction. Let us further assume that V/ is a fixed member of the Gaus¬ 
sian orthogonal (GOE) or unitary ensemble (GUE). This may be justihed 
for an environment dominated by chaotic dynamics or disorder. The Born- 
Markov approximation used above implies some coarse-graining in time, i.e. 
averaging over a small time interval. Assuming ergodicity, this averaging may 
be replaced by averaging over the random matrix ensemble. In the case of the 
Gaussian ensembles, it is then enough to use the fact that (V/jj V^ki) = ^jk 
(GUE) or (V/jj Vlki) = 5jk ^ii + ^ik Sji (GOE), to arrive at the following very 
simple master equation: 

ihdt 0e{t) = Hx Qeit) - Qeit) Hq - T { g^{t) “ ^ ^r ge{t) } . (19) 

The dissipation term of this equation had been obtained previously for a 
random matrix model where the coupling matrix is a full random matrix 


in the product Hilbert space of near and far environment [14]. Note that 
Eq. (12) implies 

f\,r{t) = trpe(i) = tTe[Mx{t) Qeit)] , Mx{t) = Uo{t)^Ux{t) . (20) 

2.3 Evolution equation for the generalized fidelity am¬ 
plitude 

To solve the master equation (19), we return to the interaction picture and 
separate off the exponential decay from the solution: 

Qe{t) = e“^* H(f) : dtVL{t) = ^ Ml{t) tTe[Mx{t) Q{t)] . (21) 

We now consider the corresponding integral equation 

n{t)=n{0) + ^[ dr M|(r) tre[MA(r)fl(r)] , ( 22 ) 

Jo 

and dehne 

0A,r(^) = tre [ Mx{t) Q{t) ] such that = e“^* 0A,r(^) • (23) 

Eventually, we are only interested in the generalized hdelity amplitude fx,r{t)- 
It is then convenient to derive an evolution equation for (j)x,r{t) directly. To 
this end, we multiply the above integral equation with Mx{t), and take the 
trace. This yields 

0A,r (t) = tre [ Ma (f) Pe (0) ] + ^ f dr tre [ Ma (f - r) ] 0A,r (r) (24) 

Jo 

The hrst term of the RHS is just the hdelity amplitude in the near environ¬ 
ment without coupling to the far environment. We will denote this function 
as fx(t) = fx,o{t)- The term N~^ trg [ Mx{t — r) ] is the same type of hdelity 
amplitude, just that here the initial state is the maximally mixed state 1 /iVe- 
That function we will denote with fx{t)- In the numerical simulations in 
Sec. 3 we will assume that Pe(0) is the same maximally mixed state, such 
that fx(t) = Thus, we may write 

4>x,r{t) = /a(^) + T [ dr fx{t - r) 4>x,r{r) . (25) 

Jo 
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Since the integral in this expression is precisely equal to the convolution of 
the functions f\(t) and 0A,r(t), the equation can be solved formally by a 
Laplace transformation [17]. It leads to 

( 26 ) 

and similarly for F{s). In this expression, the ratio F{s) [1 — L.F(s)]“^ can 
be formally expanded into a power series in L, where the inverse Laplace 
transform of powers of F{s) and F{s) yield iterated convolutions of the orig¬ 
inal hdelity amplitudes. Terminating the series at hrst order in L leads to 
the approximate expression for derived in [5]. 

In practice, it is more convenient to solve the integral equation directly 
using a numerical equal-stepsize integration scheme. The theoretical curves 
compared to random matrix calculations in the following section, are ob¬ 
tained via the trapezoidal rule. We noticed that the larger T the smaller the 
required stepsize, which does not allow a straight forward exploration of the 
large T limit. 

3 Numerical results 

In this section, we verify the integral equation (25) for the generalized hdelity 
amplitude, when the coupling to the far environment can be described by an 
unitarily invariant dissipation term, see Eq. (19). We have shown that this 
situation occurs when the dynamics in the near environment can be described 
by a random matrix ensemble, for instance in the case when the quantum 
chaos conjecture applies. 

For the numerical simulations, we use the methodology of Ref. [5], writ¬ 
ing the evolution equation for ge{t) in super-vector and super-matrix form, 
and solving the resulting system of differential equations by diagonalization. 
Since the number of equations is of order N^, we are restricted to relatively 
small dimensions, W = 50. For each case, we perform three independent 
numerical averages over nrun = 1000 realizations. This allows to estimate 
(roughly) the statistical uncertainty of the numerical results. 

To obtain a solution of the integral equation, we employ the method 
explained at the end of the previous section, which is, as far as precision and 
computer workload are concerned, equivalent to direct numerical integration. 
Hence, for not to large values of T, we can obtain highly accurate results 
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Figure 1: Generalized fidelity amplitude /A,r(t) subtracted by fx(t) for A = 
0.1 and different values of F: F = 0.01 (solid line), 0.05 (dashed line), 0.077 
(dotted line), and 0.1 (dashed-dotted line). The lines show the solutions of 
Eq.(25) and the points show the corresponding numerical simulations. 


in very short computation times. For the model considered, the hdelity 
amplitude f\{t) without far environment, would be given by the universal 
random matrix result hrst derived in Ref. [18]. However, for simplicity we 
use the exponentiated linear response formula from Ref. [19]. 

Below, we consider two different strengths of the dephasing coupling. In 
Fig. 1 this is A = 0.1 which corresponds to the cross-over regime where the 
decay of fx{t) is intermediate between exponential and Gaussian decay. In 
Fig. 2, A = 0.02 which corresponds to the perturbative regime where the 
decay of fx(t) is dominantly Gaussian. In both hgures, the time scales are 
chosen such that the Heisenberg time is at t = 27r. 

In order to compare the theoretical prediction with the numerical simu¬ 
lations on a hner scale, we will plot the difference between the generalized 
hdelity amplitude fx,rit) and the hdelity amplitude fxit) = without 

coupling to the far environment. 

In Fig. 1 we can clearly see the inhuence of the coupling to the far envi¬ 
ronment, on the coherence measured in the central system. The ehect scales 
with a = F/A, which, in this hgure ranges from a = 0.1 to a = 1. Since 
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Figure 2: Generalized fidelity amplitude fx^rit) subtracted by f\(t) for 
A = 0.02 and different values of F: F = 0.00195 (solid line) 0.01 (dashed line) 
0.02285 (dashed-dotted line) and 0.1 (dashed-dashed-dotted-dotted line). 
The dashed lines show the solutions of Eq.(25) and the points show the 
corresponding numerical simulations. 


the difference plotted is positive and increasing with a, we conhrm the effect 
discussed in Ref. [5], that increasing the coupling to the far environment, sta¬ 
bilizes the coherence in the central system. We also observe, that the theory 
obtained from Eq. (25) agrees very well with the numerical simulations for 
all values of a and for all times. 

In Fig. 2 we repeat the comparison but for the case of A = 0.02, which 
is well in the perturbative regime, where fx{t) shows a Gaussian decay. 
Again we compare numerical simulations with the theoretical prediction from 
Eq. (25), for different values of a. Here, these values range from a = 0.0975 
up to a = 5. In this range, the stabilizing effect of the coupling to the far 
environment is even stronger. However, we also observe that the theoretical 
prediction systematically overestimates the effect, the larger a the larger the 
deviation. We suspect that these deviations might be a hnite W effect. 
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4 Conclusions 


We have presented a generalized fidelity amplitude, which starts out from the 
only direct measurement scheme for the fidelity amplitude [1], which uses a 
central system as a probe for the echo dynamics in the environment. We open 
the system to be probed to an additional outer environment and consider the 
decay of coherences of the probe. This we define to be proportional to the 
generalized fidelity amplitude. While this ensures that the quantity is mea¬ 
surable, it remains to be investigated how it relates to other generalizations 
of fidelity to open systems, such as the one by Josza [20]. 

We obtain a more general master equation, Eq. (17), which may be ap¬ 
plied to integrable and non-integrable models for open quantum systems. The 
random matrix model considered here to illustrate our results, is slightly dif¬ 
ferent from the model considered in [5], since there we considered the far 
environment to be a bath of harmonic oscillators described by the Caldeira- 
Leggett model of random Brownian motion [21]. Nevertheless, its effect on 
the near environment is equivalent and can be described by the same unique 
parameter T. We derived a new exact integral equation for the generalized 
fidelity amplitude, which constitutes significant progress comparing to previ¬ 
ous work. Expanding this equation to lowest order in T leads to the approx¬ 
imate analytical formula obtained in [5]. In Sec. 3, we find that its solution 
agrees with numerical simulations over a broad range of coupling strengths. 
This confirms the general experience, that increasing the coupling between 
near and far environment protects the coherence in the central system. 
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